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I .   INTRODUCTION 


A. 


BACKGROUND 


This  paper  describes  the  concept  and  construction 
of  a  hard  wired  analog  device,  pictured  below,  which 
simulates  the  operation  of  a  gas  turbine  engine. 


Recent  work  on  gas  turbine  simulation  has  focused 
largely  on  digital  or  hybrid  computing.   To  accommodate 
the  nonlinear  elements  of  a  gas  turbine  simulation  in  this 
fashion  requires  a  perturbation  or  piecewise  linear 
technique.   A  pure  analog  approach  offers  a  continuous 
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solution  to  these  nonlinear  characteristics,  which  facili- 
tates a  simulation  in  real  time. 

At  this  time,  analog  computation  of  complex  non- 
linear systems  is  not  feasible  on  devices  such  as  the  EAI 
680  computer  currently  installed  in  the  M.I.T.  Joint 
Computer  Facility.   With  the  miniaturization  of  both  digital 
and  linear  functions  into  the  now-common  integrated  circuit 
or  IC  chip,  it  becomes  possible  to  design  extremely  compact 
analog  devices.   The  IC  component  cost  of  such  a  simulation 
is  rapidly  approaching  the  cost  of  purely  passive  electronic 
elements. 

The  modeling  approach  to  the  analog  simulation, 
physical  modeling,  follows  the  work  of  Markunas   closely. 
Physical  modeling  attempts  to  explain  the  engine  char- 
acteristics within  a  framework  of  physical  princples  and 
laws  which  describe  the  functions  occurring  in  the  real 
machine.   Simplicity  is  generally  a  strong  point  of  the 
model. 

The  model  which  has  been  implemented  in  the  analog 
simulation  represents  a  general  machine  "invented"  by 
Markunas,  who  estimated  the  necessary  characteristic 
relations  between  the  independent  variables.   Markunas 
conducted  digital  and  hybrid  simulations  of  his  "general" 
gas  turbine,  but  a  full  analog  implementation  was  not 
investigated. 
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This  paper  is  essentially  a  design  manual  which 
traces  the  development  of  a  hard-wired  analog  simulator. 
It  is  intended  to  give  the  reader  a  brief  overview  of  the 
analog  devices  currently  available,  and  their  application 
to  a  computational  circuit.   To  facilitate  future  work 
with  the  simulator,  a  detailed  user's  manual  has  been 
written. 

This  document  contains  explicit  details  relating 
to  the  construction  and  operation  of  the  simulator.   The 
user's  manual,  the  simulator,  and  associated  hardware 
are  under  the  supervisor  of  Professor  Henry  M.  Paynter, 
Department  of  Mechanical  Engineering. 

B.      DESIGN  CONCEPT 

The  following  characteristics  were  desired  in  the 
hard-wired  simulator: 

1.  Real  time,  non- linear  third  order  simulation 
of  the  dynamics  of  a  single  shaft  gas 
turbine  engine,  from  start-up  to  overload 
conditions. 

2.  Cost  to  be  under  $600. 

3.  Size  of  the  device  such  as  to  be  readily 
portable  and  convenient  to  set-up  and  operate. 

4.  Layout  and  operation  of  the  simulator 

to  be  comparable  with  typical  industrial 
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and  military  engine  control  consoles. 

5.  Accuracy  of  the  steady-state  simulation  to 
be  considered  of  secondary  importance 
compared  with  the  representation  of  large 
scale  transient  effects. 

6.  The  ability  to  investigate  various  fuel 
control  schemes  in  a  closed  loop  control, 
and  the  ability  to  operate  with  different 
torque-speed  load  characteristics. 

7.  Utilization  of  the  general  purpose,  low 
cost  IC  operational  amplifier  in  a  hard 
wired  configuration  to  accomplish  all 
mathematical  operations. 
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II.       MODELING 

A.  GENERAL  APPROACH 

The  gas  turbine  model  selected  by  Markunas  will  be 
briefly  explained.   The  modeling  approach  is  broken  down 
into:   steady  state  and  dimensional  considerations, 
frequency  analysis,  and  dynamic  equations. 

Markunas1  estimated  performance  parameters  and 
the  dynamic  equations  are  combined  in  state-variable 
format.   These  equations  appear  in  Appendix  II  and  form 
the  basis  for  the  analog  signal  flow  used  in  the  design. 

B.  STEADY  STATE  AND  DIMENSIONAL  ANALYSIS 

3 
Horlock  establishes  non-dimensional  functions  in 

describing  the  performance  of  compressors  and  turbines. 

A  set  of  relations  are  developed  to  describe  the 

performance  of  a  turbomachine  in  its  simplest  terms. 

In  the  subsequent  development,  gas  turbine  state 

points  are  used,  numbered  in  the  following  convention,  with 

total  values  indicated  by  a  zero  in  front  of  the  state 

point: 
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£>o-v/\<j^r 


The  following  symbols  are  used  for  dimensional 
analysis  in  the  compressor: 


Symbol 

Ryz 

T0-2 

M 


Compressor 

Inlet  stagnation  pressure 

Inlet  stagnation  density 

Outlet  stagnation  pressure 

Inlet  stagnation  temperature 

Efficiency 

A  characteristic  linear  dimension 

Stagnation  temperature  rise  through 
compressor 

Rotational  speed 

Mass  flow  rate  of  gas 


Dimen 

sions 

M/ 

„T^ 

n/ 

u^ 

M  / 

*LT* 

<3 

hl/T 
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Symbol               Compressor  Dimensions 

S\L  Molecular  wt  of  gas 

^<K  Absolute  viscosity  of  gas  r^/iT 

^P  Specific  heat  at  constant  pressure  U  /T  ^ 

*y  Ratio  of  specific  heats  — 

£*s  Gas  constant  U  /T  O 

^G'Z  Speed  of  sound  of  entering  gas  L./«" 

\jLm  Torque  input  M  L.  /T" 

In  the  compressor,  if  fca*&  '^   ix   an^  /\T©  are  chosen  as 
variables  of  interest,  one  can  write 

Assuming  a  perfect  gas,  fixing  the  gas  constant ,  CL,  choosing 
"^  /  an<^  Tc^Z  thus  constrains  Qc^.  nCj-z,  can  be  expressed  in 
terms  of  (<  ,  fc^z.  IcvS   ;  the  original  relation  becomes 

Taking  non-dimensional  groups, 

With  constant  values,  V<  and  t<f  can  be  eliminated: 

Multiplying  the  last  two  non-dimensional  parameters  to 
create  a  new  parameter, 
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It  can  be  seen  that  the  third  parameter,  ^TTS   is  a 
"flow  Reynolds  number". 

Assuming  that  the  compressor  operates  above  some 
critical  Reynolds  number,  and  that  the  dimensions  of  a 
particular  machine  are  fixed,  the  equations  reduce  to: 


s"\>%.    ,^k  =  W-£L  > 


In  order  to  present  meaningful  performance  relations, 
compressor  data  should  be  corrected  following  the  above 

parameters. 

4  . 
Identically  for  a  turbine,   it  can  be  shown  that 

the  same  relations  hold  and  that  the  parameters  of 

interest  in  turbine  performance   are  {V\  y  s^l   .  1  "Z.    )&>\&- 

All  inlet  conditions  assume  values  at  state  point  4,  and 

all  exit  values  are  at  state  point  5.   Torque  across  the 

turbine  is  i*z  .   As  shown  before, 


The  characteristics  for  a  given  turbomachine  can 
now  be  specified  as  functions  of  dimensionless  or  pseudo- 
dimensionless  groups: 
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Compressor: 


Turbine : 


Po4 


£     (Ik*        *       \ 


p«z 

Mj^-2    / 

Fc*  ~  +(     pc*     >  'iTVA 
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C.  FREQUENCY  ANALYSIS 

Modeling  of  the  processes  which  take  place  in  a  gas 
turbine  can  generally  be  classified  as  mechanical,  fluid, 
or  thermal  in  nature.   To  reduce  the  scope  of  the  problem 
to  a  more  manageable  field,  only  the  most  significant 
energy  domains  and  dynamics  were  considered.   Markunas 
defines  three  dynamic  regimes  for  the  gas  turbine: 

1.  Low  frequency,  (f  -s^  .1  hz)  where  rotary 
inertias  and  thermal  capacitances  dominate 
the  dynamics. 

2.  Medium  frequency  (f ^  1  hz)  where  rotary 
inertias  and  fluid  capacitances  dominate. 

3.  High  frequency  (f^t"  100  hz)  where  the  fluid 
dynamics  are  dominant. 

Markunas  established  that  rotor  dynamics  and  fluid 
capacitances  were  the  primary  influences  on  gas  turbine 
dynamics.   Thermal  capacitances  and  fluid  inertias  were 
neglected. 

D.  DYNAMIC  EQUATIONS 
1.   Assumptions 

The  gas  turbine  model  is  assumed  to  operate  with  a 
perfect  gas  in  one-dimensional  flow.   Components  are 
treated  in  lumped  parameter  fashion. 


-17- 


2.   Rotor  Inertia 

For  the  rotor,  a  summation  of  torques  give  the 
following: 


tt  T4* 

where 

«k 

S^Cr       =   angular  acceleration, rpm/sec 

2 
XLc^,       =  rotary  inertia,  f  t-lbf-sec 

"\"       =   load,  compressor,  or  turbine  torque, ft-lbf 

3.  Fluid  Capacitance 

For  the  one-dimensional  lumped  parameter  fluid,  the 

following  symbology  is  used: 

£>  =   fluid  density,  lbm/ft3 

y{  =   axial  distance  along  duct,  ft 

S/y       =  mean  axial  velocity  of  fluid,  ft/sec 

*£$        =   ratio  of  specific  heats  =   Cp/Cv 

Cp       =  specific  heat  at  constant  pressure 

Cv/       =  specific  heat  at  constant  volume 

\**\  =  Mach  number  =   V*  /  ^< 

Cv       =  Local  sonic  velocity  =  >/  J*-  ^  ^ 
<J        =   internal  energy  of  the  fluid  in  ft-lbf/lbm 
^>       =  enthalpy  of  the  fluid  in  ft-lbf/lbm 
V^L     =   lumped  fluid  volume 
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A  general  mass  conservation  equation  for  the  one 
dimensional  fluid  stream  may  be  written: 


Integrating  over  a  constant  area  duct,  with  fluid  pro- 
perties constant  at  a  given  axial  distance  or  dimension 


-^ — ■   «E~  Ci       by  definition. 
Integrating  over  the  length  of  the  duct 


^5r 


JfMy      *    ^^Cf^V^cU  so 


^.^.V^    *.    ^/W* 


^f  A,ds 


is  approximated  by  an  average  density, 
f   •     Lumped  fluid  volume  (V^L.)  since 
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where  the  subscipts  1  and  2  identify  inlet  and  outlet 
conditions  across  the  engine,  respectively. 

Writing  a  similar  equation  for  the  conservation  of 
momentum  neglecting  viscous  shear  and  body  forces, 


-£*+- 


Integrating  over  area, 


^s 


Integrating  over  length  of  the  duct 


4  (4-0  - 


*<*)- 


<5r 

VyAv 


r—    -a 


4-   PA, 


PN~1 


_   Q 


-=-  o 


I 


L- 


*3<-  J-z- 
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^5T 


x^)=[^Xp^*^^J 


H-~ 


Since   total  pressure 


.   *,  p^^^rf)— 


Assuming  V  =  1.4,  and  M  £  .5,  the  momentum  equation  can 
be  written: 


4-X^)~   *£fa.-*^ 


Similarly  for  energy  conservation  with  no  external 
work  or  heat   addition, 


4-    ^L. 


Integrating  over  area  and  length, 


U 


=.  o 
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-L 


( voo  4 


3t- 


W 


VT\  YAqs 


where 


Since 


K.= 


+    ""2#C- 


Substituting , 


"2. 


■*-  ^Z  •=.    cvT  (j_  *  *(V~i)  ^  ^ 


-2<jc 


'£. 


where 


- — =■  \U  «V-4M 

With   b'A  ~  -  ^  the  result  can  be  simplified  to: 


E.      STATE  EQUATIONS 

The  three  conservation  equations  are  now: 
1.   Mass 


iC^-  ^slC^-^ 


HF^   ; 


VOL. 
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2 .      Momentum 


^(^  -  ^  (?^~  &?) 


3 .      Energy 


-gp  CfO  *  ^(^xTol  -  ^ZT< 


o-tf 


Using  the  definition  of  a  polytropic  process, 

_  n 

p\/  -=  cot^T",  it  follows  that 

rs 


In  order  to  get  the  mass  and  energy  equations  into  more 
convenient  form, 


^•)  -fetes)  4*  ^ 
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With  the  result  that  the  equations  now  be  more: 


4t  c&  =  *&(**-  *J 


"3f 


Where  1  and  2  represent  inlet  and  outlet  across  the 
duct.  Because  of  lumped  parameter  considerations,  machine 
model  parameters  can  be  substituted  as  follows : 


One  Dimensional 

Machine  Model 

<v\* 

a\z 

W^ 

~r*> 

Tc^ 

15*. 

To-s 

TS'Z. 

Tc*4 

P« 

Po* 

Substituting  notation  for  the  first  and  third 
differential  equations, 
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T^   =  (f^^^^T^^^V^T^ 


Introducing  the  torque  equation  as  the  third  system 
equation , 


^fr*  -it  -fCI 


where 

^LUV       =   fuel  lower  heating  value,  Btu/lbm 

W\p        =  fuel  flow,  lbm/sec 

'Tt- 
*  ^        =   load  torque,  ft-lbf 

r^a^   TW^   <VhO  t*N^-    are  chosen  as  variables  of 
convenience. 

For  an  isentropic  process,  n  =  1.4.   Markunas  felt 
that  n-l/n  should  be  replaced  by  1.0  in  the  temperature 
equation. 

F.      ENGINE  DATA,  ASSUMPTIONS,  AND  CHARACTERISTICS 

1.   Introduction 

In  Markunas'  analysis  a  single  shaft  machine  was 
invented  from  a  number  of  sources.   Several  real  world 
machines  exist  which  have  design  parameters  close  to 
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Markunas's  engine.   The  recent  Garrett  GTPF  990  engine  in 
its  locked  shaft  mode  is  a  reasonable  example. 

2.  Design  Point  Specifications 

Pressure  Ratio  5.0 

Compressor  Isentropic  Efficiency  1280°F 

Turbine  Isentropic  Efficiency  88% 

Compressor  Discharge  Mass  Flow        100  lbm/sec 
Inlet  Pressure  14.175  psia 

Inlet  Temperature  80  F 

Rotor  Speed  7200  rpm 

Output  6948  HP 

Markunas  plotted  characteristic  curves  for  his 
compressor  and  turbine,  which  appear  in  Appendix  I. 
Markunas  then  mathematically  "fit"  each  characteristic 
curve.   These  approximations  follow  the  functional 
relations  established  for  non-dimensional  and  pseudo 
dimensionless  groups. 

3.  Compressor  and  Turbine  Characteristics 
Markunas  defined  the  following  variables  and  char- 
acteristics for  his  compressor: 

PRC   ~      Fo-i/fe-Z.  pressure  ratio 
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TRC  -=        Xo^  /r< 


0-2  temperature  ratio 


WAw/r  -  ^^^^"^//V^^sili-zX         scaled  compressor 
^Wk-      fU.  A  "p™  7        mass  flow 

XNC   -   N<ar,_  ,/f4<sr  \  scaled  compressor 


EC      *■         x.o     4-    XKlC   *■  -S^CxMC^) 


PC     -  4--0        l(r        ppC  4:  gc 

-z-o    ip     p££  sec 
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QrEC     -  ^^r 

cc 

L  BC 

I—  F° 

EFFC     =       DC  i  -  (kc  -  PCC\ 

L       NEC-1     /. 

(PRC  )     -  L 


Similarly,  for  the  turbine 


P^T       -    Po<r/^ 


*"  4$L) 


A»  I        -     1  .  OOQ^ 


gr  ■*     -^a.   -*- 
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FT 
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MM 
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Additional  Machine  Assumptions 

1.     R^S  =  f^2  -  F^X      (no  duct  loss) 

,~  ^  (no  combustor  losses) 

Pcrs  -  P*-* 

The  machine  characteristics  are  combined  with  the 

differential  equations  to  produce  a  set  of  system  state 

equations.   The  expanded  equations  appear  in  Appendix  II. 


A. 
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III.   THE  OPERATIONAL  AMPLIFIER 


INTRODUCTION 


The  fundamentals  of  the  operational  amplifier,  (OA) 
are  briefly  presented  to  establish  the  necessary  analog 
"tools"  to  proceed  to  a  simulator  design. 


B 


BASIC  OPERATION 


An  ideal  OP  amp  can  be  modeled  by  the  following 


with  no  differential  voltage  between  positive W  )     and 
negative  V V  /inputs,  both  emitters  are  equally  biased,  and 
collector  currents  Qi-  and  C^2are  equal.   The  differential 
voltage  output  between  1  and  2  is  also  zero.   Driving  V 
more  positive  with  a  signal  will  cause  the  current  through 
Q^to  increase.   The  current  through  Q^  must  decrease  by  a 
proportional  amount,  so  that  in  the  differential  voltage 
output  between  1  and  2 ,  1  gets  driven  more  negative  with 
respect  to  its  initial  voltage,  and  2  becomes  more  positive. 
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Thus  for  a  small  signal  applied  to  one  of  the  input  terminals, 
a  high  voltage  gain  and  current  "buffering"  takes  place  as 
a  result  of  the  transistor  characteristics. 

With  the  ideal  OA,   the  input  current  is  essentially 
zero,  and  the  voltage  gain  is  assumed  infinite. 

C.      THE  INVERTING  OA 

Looking  at  the  following  diagram, 

pA/WV^Tb 


Xo," 


±   Jp-4*) 


It  is  assumed  that  the  impedance  of  the  amplifier 
"He-3**0  ,  therefore  Tc%0_ 

Closed  Loop  Gain,  from  a  summation  of  currents: 


Xb       a      (Vo-vO/s- 


Cvr  -  vO/ 


.l 


r 


<X 


1       * 


/^    ♦ZaS] 


^Zi.//^  +z 
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If  a,  the  amplifier  gain,  is  large, 


X. 


It  can  therefore  be  seen  that  the  inverting  OA  inverts 
and  provides  a  gain  proportional  to  the  input  and  the 
signal  feedback  resistances. 


D. 


INVERTING  SUMMER 


By  joining  multiple  input  currents  to  the  input  side 
of  the  OA,  an  inverting  summer  can  be  constructed  with 
variable  gains  on  each  input. 


r-AAMA-i 


s/o 


By  current  summation, 


Vc  /h- 


Vl/^j,.   -  V-Z./W 


■;£ 


v. 


"3r. 


O 


-  S/- 
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E. 


INTEGRATOR 


Adding  a  capacitor  to  the  feed  back  loop  of  the  OA 
adds  the  following  constitutive  relation  to  the  feedback 
loop: 


<-\$- 


Looking  at  a  typical  circuit, 


tt 
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x^  ~~  -r< 
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F. 


LOGARITHMIC  GENERATOR 


Using  a  diode  in  the  feedback  loop  enables  the 
logarithmic  nature  of  the  diode  voltage  -  current  relation 
to  be  included  in  the  loop.   The  diode  characteristic 
can  be  modeled  by  the  following  equation: 


v 


L 


IN 


1 


XT* 


-   ^/o  Lu  (x*   v  E  ") 


where : 


"Ho 

T 

V 


=  junction  saturation  current 

=  junction  saturation  voltage 

=  bias  current 

=  bias  voltage 


A  typical  characteristic  appears  as: 


iv\A» 


To 
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It  is  important  to  note  that  the  saturation  voltage  ^o  is 
directly  proportional  to  the  absolute  temperature. 

A  typical  log  circuit  can  be  constructed  as  follows 


V, 


A/WW 


A  diode  -  connected  transistor  is  used  to  perform 
the  diode  function,  as  transistor  emitter-base  junctions 
have  more  dynamic  current  range  than  diodes  in  general. 
By  the  following  diagram,  an  input  \V  across  £1.  causes 
a  current  flow  through  the  "transdiode" .   The  bias  voltage, 
appearing  as  ^O)  varies  as  the  log  of  V»"  . 


Coe&€"»^T 


S/oun^  fe 
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G.      ANTILOG  GENERATOR 

In  a  similar  fashion,  input  voltages  or  currents  may 
be  exponentiated  in  the  following  circuit: 


V/i' 


T^ 


— www- 


«■ 


-o 
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IV.   ANALOG  COMPUTATION 

A.  INTRODUCTION 

To  implement  the  non-linear  system  equations  the 
following  functional  operations  are  performed:   addition, 
subtraction, integration,  multiplication,  division, 
exponentiation  to  a  fixed  power,   and  exponentiation  to 
a  variable  power.   The  first  three  analog  operations  have 
been  briefly  discussed.   To  implement  the  higher  functions, 
the  use  of  logarithmic  circuits  must  be  made. 

B.  THE  PRACTICAL  LOG  CIRCUIT 

There  are  at  least  three  major  considerations  in 
the  construction  of  a  practical  log  or  antilog  circuit: 
cost,  complexity,  and  accuracy.   Three  different 
architectures  are  considered  to  accomplish  the  log  function. 

The  first  choice  is  the  previously  discussed  simple 
log/antilog  circuit.   Using  low  quality  components,  cost 
of  the  log  circuit  is  under  one  dollar.   Thermal  stability 
is  considered  poor.   A  test  circuit  and  transdiode 
characteristics  are  shown  in  Appendix  III.   Test  data 
indicated  that  the  accuracy  and  thermal  properties  of 
this  circuit  were  marginal. 

A  second  choice  is  a  compensated  log  circuit  as 
diagrammed  below,  which  offers  increased  thermal  stability 
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with  the  use  of  a  thermistor: 


14 


Voetp 


Q: 


-^wwv-^ 


r©    r 


a/VWH 


fUg^S**^. 


In  the  diagram,  the  collector  of  $z  is  kept  at  virtual 
ground  by  the  OP  amp,  and  the  collector  current  through 
C^z  is  equal  to  a  constant  reference  current  provided 
by  vfc-e^  .   It  can  be  shown  that  the  output  voltage  is  of 
the  form   Vo  —  ~*  L-C<**1«>  v--,—  '  )    where  the  thermistor 
provides  temperature  compensation.   One  of  the  drawbacks 
to  this  scheme  is  that  the  number  of  components  necessary 
to  generate  the  log  function  has  gone  from  3  to  at  least 
9,  with  proportional  space  and  cost  increases.   This 
configuration  was  not  evaluated. 

The  third  and  most  expensive  choice  is  to  utilize 
a  dedicated  logarithmic  device.   The  particular  unit 
tested  was  a  Texas  Instrument  log  amplifier.   This  is  a 
dual  log  generator  on  a  single  DIP  package,  at  a  cost  of 
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approximately  $2.50.   With  appropriate  circuit  manipulation, 
the  desired  log/anti-log  operations  can  be  performed.   The 
TI  log  amplifier  was  evaluated  in  depth.   Results  are 
documented  in  Appendix  III.   The  log  amplifier  was  the 
foundation  of  the  simulator  design.   Other  dedicated  log/ 
antilog  modules  exist  which  offer  increased  dynamic  range, 
temperature  stability,  and  accuracy,  but  alternative 
devices  were  an  order  of  magnitude  more  costly  than  the 
TI  log  amplifier. 


C. 


EXPONENTIATION  -  ALL  CASES 


In  order  to  achieve  either  fixed  or  variable  power 
exponentiation,  the  use  of  logarithms  is  basic.   The 
signal  flow  for  a  fixed  power  exponent  (power  function 
generator  or  PFG)  circuit  is  diagrammed  below: 


yvww 


LjCBj 


(Z^ 


-WWW  -  - 


*  i. 


--A/VWV 


By  a  judicious  selection  of  values  for  \£L  and  R^ 
current  can  be  multiplied  or  divided  using  an  operational 
amplifier  in  the  ratio: 


-39- 


l 


To  be  able  to  exponentiate  to  a  variable  power 
requires  a  cascaded  log/antilog  arrangement.   The  use  of 
a  segment  of  the  state  equations  is  made  to  demonstrate 
the  concept.   It  is  desired  to  generate  the  following: 

where  K.  represents  the  scaling  of  several  fixed  parameters, 
and  X^^PRX  }   8^   and  ^  are  all  problem  variables. 
The  circuit  selected  must  process: 


wc 


As  a  result  of  thermal  instabilities,  the  interior  section 
of  the  signal  flow  was  replaced  with  a  four-quadrant 
multiplier. 

The  multiplier  essentially  compresses  the  separate 
functions  into  a  single  chip.   By  circuit  manipulation, 
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a  multiplier  can  also  divide,  square,  and  take  square 
roots.   A  range  of  multiplier  devices  currently  exist, 
with  the  price  of  the  least  expensive  about  $3. 

Test  data  on  the  multiplier  used  in  the  simulator 
was  collected  in  Appendix  III.   To  allow  as  compact  a 
circuit  as  possible,  the  multiplier  was  frequently 
utilized  in  its  dividing  and  exponentiating  modes. 
Complete  device  specifications  are  included  in 
Appendix  VI. 
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V.       PRELIMINARY    DESIGN 

A.  INTRODUCTION 

The  simulator  design  proceeded  in  a  parallel  fashion. 
In  an  iterative  process,  estimations  were  generated  for  the 
electronic  layout,  which  allowed  the  cabinet  and  control 
panel  to  be  sized  and  designed. 

The  simulator  is  configured  to  resemble  a  real-world 
engine  control  panel.   Additional  concepts  for  the  per- 
ipheral design  which  were  not  incorporated  were: 

1.  Proportional  fuel  control  schemes  with 
variable  gains. 

2.  Torque-speed  loads  characteristic  of 
pumps,  propellers,  generators,  and  other 
real  world  devices. 

3.  Overspeed  alarms  and  stall  audio  noise 
circuits. 

B.  ELECTRONICS 

1.   Logic  Circuit 

The  logic  processing  of  the  electronic  circuitry 
was  broken  down  into  3  separate  grid  structures,  which 
appear  as  design  drawings  in  Appendix  V.   Use  was  made 
of  bread-board  panels  to  avoid  the  inflexibility  of  other 
techniques.   Concurrent  with  testing  various  mathematical 
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operations,  the  circuit  design  was  modified  to  balance  the 
competing  parameters  of  simplicity,  size,  accuracy,  and 
cost.   The  appropriate  equation  variables  were  normalized 
about  their  design  point  values,  and  then  scaled  to 
provide  a  circuit  signal  that  was  approximately  75%  of 
the  maximum  rates  input  for  the  device.   The  steady 
state  design  point  values  for  all  variables  were  established 
by  a  DYSYS  program  and  used  for  the  normalization. 
Additional  notes  on  the  wiring  techniques  and  breadboard 
layout  are  included  in  Appendix  VI. 

2.   DYSYS  Simulation  and  Scaling 

A  DYSYS  simulation  of  the  system  equations  was  run 
as  it  appears  in  Appendix  IV.   Significant  comments  on 
the  program  are : 

1.  The  conditional  value  of  FC  was  changed 
to  a  fixed  value  of  3.0  to  facilitate  the 
circuit  logic. 

2.  Markunas  estimated  a  parabolic  "pump-type" 
load  torque,  which  appears  as  .0000977  (3) **2. 

3.  The  modeling  parameter  DT  is  called  TT  in 
the  program  to  avoid  confusion  with  the 
DYSYS  time-step  variables,  DT. 
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The  scaling  procedure  started  with  the  following 
listings : 

a.  Constants 

Tb-z.  STi-V  (Z 

b.  Machine  Characteristics 


(fe) 


^l^ 


(m-^sj  i"*-  ) 


_L.i~sa 


f^-2   OfeS.U«*4 


^oS^ 


'Oe:»iCr«a 
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c.  Input  Variables  at  Design  Point 

d.  State  Variables  at  Design  Point 

Tc*4         «s  lT7'^r   6L 


The  DYSYS  simulation  was  run  to  verify  Markunas's 
simulations.   Use  was  made  of  the  program  to  assist  in 
adding  voltage  and  current  values  to  the  circuit  design, 
and  to  evaluate  the  effect  of  fixing  the  discontinuous 
constant,  FC.  FC  was  conditionally  assigned  the  value 
of  2  or  4  in  Markunas  work.   The  steady  state  effect  of 
assigning  FC  as  3  was  investigated  and  adopted. 

The  DYSYS  simulation  was  also  run  off  -  design 
to  generate  a  test  case  for  the  logic  circuit.   The 
simulation  was  run  at  approximately  50%  fuel  flow  and 
steady-state  values  were  generated  for  all  elements  of 
the  equations. 
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Some  difficulty  was  encountered  in  determining  the 
correct  value  for  the  lumped  parameter  fluid  volume.  The 
value  which  appeared  to  be  consistent  with  Markunas1  data 

was  144  ft  .   An  approximation  based  on  real  world  data 

3 

gave  15  ft  as  a  more  likely  figure.   Since  the  volume 

represents  a  time  constant  scaling  for  the  fluid  capacitance, 
it  has  no   significance  for  the  DYSYS  steady  state  values. 
Changes  in  the  volume  assumed  can  be  readily  investigated 
in  the  simulator  dynamic  response. 

The  scaled  voltage  values  and  the  system  variables 
they  represent  are  drawn  on  the  3  circuit  plans  in 
Appendix  V. 

C.      PERIPHERALS 

The  main  elements  of  the  remaining  design  work 
included  the  chasis,  power  supply  and  distribution, 
instrumentation  and  control  hardware,  and  breadboard 
architecture.   The  design  drawings  and  sketches  included 
in  Appendix  V  illustrate  the  final  peripheral  design. 
The  photograph  below  shows  the  internal  layout  of  the 
simulator. 
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Appendix  VI,  specifications,  contains  additional 
data  on  the  peripheral  equipment,  its  wiring,  and  operating 
parameters. 
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VI.       IMPLEMENTATION   AND    CHECKOUT 

A.  PROCEDURE 

To  establish  the  characteristics  of  the  various 
devices  and  to  investigate  the  methodology  of  the  design, 
10  tests  were  conducted.   The  significant  results  are 
discussed  here  in  the  context  of  the  design  process.   All 
test  data  is  included  in  Appendix  III. 

B.  LOG  GENERATION 

Tests  1  and  2  determined  the  voltage  -  current 
relations  for  the  transdiode  -  connected  OP  amp,  using 
both  2N2222  and  2N2904  transistors.   It  became  apparent 
that  the  allowable  voltage  values  of  the  log  circuit 
were  so  narrow  that  only  current  variations  could  be  used 
for  the  signal,  over  about  a  1  1/2  decade  range  of  values. 

C.  POWER  EXPONENTIATION 

Tests  3,  4,  5,  and  6  were  designed  to  varify  the 
fixed  power  exponentiation  circuitry.   A  workable 
configuration  was  wired  in  test  6,  using  a  current  - 
dividing  resistive  network. 

Tests  7,  8,  9,  and  10  were  conducted  to  validate 
the  variable  -  power  exponentiation  procedure.   The  inability 
of  the  transdiode  -  connected  OP  amp  to  provide  a  stable 
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cascaded  log  circuit  led  to  the  investigation  of  a 
4  quadrant  multiplier  in  Test  8. 

The  Exar  2208  multiplier  characteristics  were 
investigated  to  generate  hookup  diagrams,  trimming 
values,  and  allowable  voltage  values. 

With  Test  9,  a  simple  log/multiplier/simple 
antilog  circuit  was  wired  and  observed.   Stability  was 
significantly  improved,  but  was  still  considered  below 
the  accuracy  level  desired. 

In  Test  10,  a  Texas  Instrument  441  log  amplifier 
was  tested  to  determine  its  log/antilog  characteristics. 
The  log- amp  was  incorporated  into  the  variable  expon- 
entiation circuit  and  the  stability  evaluated.   The 
steady  state  accuracy  of  the  log  amplifier  in  the 
variable  exponentiation  circuit  was  somewhat  better  than 
the  simple  log  configuration.   Since  the  dedicated 
log  amplifier  offered  better  compensation  characteristics 
than  the  simple  log  configuration,  it  was  substituted 
at  all  points  in  the  circuit  design. 

D.      GENERAL  CIRCUIT  DESIGN 

Concurrent  with  the  test  sequence,  the  overall 
circuit  design  was  changing  as  new  test  data  was  generated. 
The  more  significant  problems  which  caused  iteration  of 
the  circuit  are  briefly  discussed.   One  immediate  limitation 
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exists  using  the  TI  log  amplifier  module.   Looking  at 
"TL  441  Transfer  Characteristics"  included  in  Appendix  III, 
the  log  amplifier  only  has  a  useable  output  range  of  0  - 
450  millivolts.   This  can  be  doubled  by  wiring  an  OP  amp 
to  the  input  side  of  the  log  amplifier  to  use  the  dual 
input  capability  of  the  device.   Still,  the  output  accuracy 
is  limited,  since  various  offsets  for  multipliers  and  OP 
amps  can  easily  reach  10  millivolts. 

As  a  general  rule,  it  was  decided  to  trin  the 
circuit  at  a  limited  number  of  points,  rather  than  dealing 
with  the  offsets  required  for  each  component. 

It  was  intended  in  the  final  circuit  design  to 
operate  all  components  at  one  dual-polarity  supply  voltage. 
This  choice  nominally  would  be  +  8  volts,  the  maximum 
rating  of  the  most  sensitive  device,  the  TI  log  amp.   Test 
10  demonstrated  that  significant  nonlinearities  were  intro- 
duced into  the  multiplier  characteristics  if  the  supply  voltage 
were  reduced.   Three  solutions  were  postulated: 

1.  Run  the  multipliers  at  +  15  volts,  and  the 
OP  and  log  amps  at  +  8  volts. 

2 .  Compensate  the  multiplier  to  have  linear 
characteristics  at  +  8  volts  supply. 

3.  Replace  the  multiplier  modules  with  log 
amplifiers  to  process  the  same  function. 
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The  first  choice  was  selected  in  an  effort  to  keep 
the  circuit  design  as  simple  as  possible,  in  as  much  as 
the  single  -  voltage  supply  was  more  of  a  packaging 
consideration  than  an  electronic  one. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  FEASIBILITY  OF  THE  DESIGN 

From  the  test  data  collected,  a  real  time  simulation 
can  be  conducted  with  the  logic  circuits  wired  as  designed. 
The  techniques  developed  for  the  implementation  of  non- 
linear differential  equations  should  be  applicable  to  a 
wide  range  of  system  models. 

Markunas  defined  his  engine  in  essentially  non 
dimensional  or  pseudo-dimensionless  groups.   It  was  not 
possible  to  investigate  the  ability  of  the  simulator  to 
represent  different  geometric  machines,  however,  it  is 
felt  that  future  work  will  verify  the  utility  of  the 
simulator  and  the  nonlinear  model. 

B.  DESIGN  REFINEMENTS 

The  simulator  as  designed  could  be  significantly 
enhanced  with  more  general   torque  and  fuel  control 
schemes.   If  the  simulator  could  demonstate  sufficient 
generality  in  modeling  the  behavior  of  real  world  gas 
turbines,  an  autonomous   "observer"  circuit  could  be 
configured  for  a  wide  range  of  engines.   This  observer 
device,  using  the  logic  developed  in  the  simulator,  could 
be  an  extremely  compact  unit.   It  could  provide  improved 
control  of  fuel,  speed,  temperature,  and  pressure 
dynamics.   Such  a  circuit  could  be  incorporated  in  one 


-52- 


or  two  monolithic  LSI  chips,  to  be  externally  trimmed  for 
a  particular  engine  or  operating  condition. 
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MACHINE  CHARACTERISTICS 
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APPENDIX   II 


STATE  EQUATIONS 
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APPENDIX  III 


Test  Data 


A. 


TEST  1  AND  2 


19 


From  reference  data,    the  transdiode  voltage  current 
characteristics  was  estimated  for  various  values  of 
saturation  current.   The  log  circuits  diagrammed  subsequently 
for  the  2N2222  and  2N2904  transistors  demonstrated  low 
base-emitter  saturation  currents,  thus  making  the  exponential 
characteristic  curve  more  abrupt. 
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Test  Data 

1  =  *. 

(EQ/26Mv) 
e 

EQ(Mv) 

I  (MA) 
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8.87  x  10 
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Test  Circuit 
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TEST    2 


CIRCUIT   DIAGRAMS   AND   DATA 
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B. 


TESTS  3-6 


The  fixed  exponentiation  of  a  variable  was  evaluated 
in  Tests  3-6.   Tests  3,  4,  and  5  were  failures  as  the 
essential  concept  of  current  division  to  achieve  the 
exponentiation  was  not  understood.   Test  6  demonstrates  the 
use  of  a  resistive  network  which  accomplished  the  current 
manipulation.   It  is  now  apparent  that  as  a  result  of  the 
characteristics  of  the  simple  log  configuration,  current 
division  can  be  easily  accomplished  by  assigning  the 
appropriate  resistor  value  to  the  antilog  input.   This 
is  diagrammed: 
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TEST  3 


CIRCUIT  DIAGRAMS  AND  DATA 
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TEST  4 


CIRCUIT  DIAGRAMS  AND  DATA 
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TEST    5 


CIRCUIT    DIAGRAMS    AND   DATA 
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TEST  6 


CIRCUIT  DIAGRAMS  AND  DATA 
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C.      TEST  7 

It  is  desired  to  validate  the  circuit  segment  for  the 
operation  P£ 


3C 


The  operational  sequence  is  as  follows: 


$■-*& 


'< 


The  circuit  described  in  the  test  data  was  set  up  to 
evaluate  the  feasibility  of  the  operation. 
The  desired  operating  characteristic  of  the  cascaded 
logger  -  antilogger  is  picutred  as  follows: 
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3^-^ 
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The  first  log  geneator  establishes  a  voltage  - 
current  relation  based  on  the  diode  and  2 OK  resistor. 
Consider  the  second  log  generator  to  be  superimposed  in 
the  first  with  a  3K  resistor  and  a  second  diode. 

The  cascaded  antilog  -  configured  amplifiers 
displayed  gross  thermal  sensitivity.   The  simple  OP  amp  - 
transdiode  configuration  is  unsuitable  for  this  cascaded 
operation  without  additional  temerature  compensation. 
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TEST  7 


CIRCUIT  DIAGRAM 
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D. 


TEST  8 


i.   Purpose 

Test  8  was  conducted  to  verify  the  variable  exponent 
part  of  the  analog  circuitry,  using  the  four  quadrant 
multiplier. 

The  equation  of  interest  is  written  below: 


LC^Qr- 


MOi_A\fU^<2. 


/V^ry-LOCi- 


-> 


M- 


ii.  Methodology 

The  functional  method  used  to  exponentiate  to  a 
variable  power  was  to  regulate  the  multiplier  output 
current  in  the  ratio  called  for  by  the  exponentiation.   For 
example : 


IF   CC 


i_  ^  ex 


.  oih 
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For  multiplier  scaling  the  integral  buffer  OP  amp 
in  the  2208  was  used  for  the  current  division  with  a 
resistive  bridge  as  follows, 


Muir>£Ue<2  i-  T^ov^-fciZ  o^ 


r-A/WW-i 


such  that  i-,  =  ±2    and  the  output  voltage  of  the  multiplier 
is  thus  controlling  the  current  input  to  the  antilog 
amplifier  transdiode. 

The  multiplier  output  is  nominally: 


i-O 


V- 


4l  ^.^^   v~yW 


&j£S6»;  (s«^c 


It  was  determined  that  a  resistive  bridge  must  be 
used  for  the  multiplier  input  from  logging  OA,  sketched 
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as   follows: 


^vwvL  - 


|MOL.v»fLvt<Z 


-MWh- 


n\^ 


cc 


so  by  modulation  of  CC  input,  the  desired  current  ratios 
across  multiplier  can  be  achieved. 

iii.  Data 

A  5  minute  observation  of  the  steady  state  output 
voltage  was  made  to  investigate  the  circuit  thermal 
stability. 
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10  Second  Samples  -Volts  x  10 


-<£>"3>o 

.  fc>*3^ 

■    <5^>"S 

.^4^ 

.  fo"2ST 

.  <^rr  1 

-  fc€T7 

.  <s^ 

.  ^44 

.  t»si 

.  'cTS'S 

.  ^-cx 

-  <b>^0 

,  &S  \ 

.•€^4. 

.G**\l 

.<^L 

.  <b*S"Z 

26   Samples 

-84- 


Ex&£.   Z2Q&    MumpuER 
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E. 


TEST  9 


i.   Purpose 

To  investigate  the  Texas  Instrument  TL  441  log 
amplifier. 


ii.  Log  Configuration 
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The  output  characteristic   V^y  ss  2  L-oC*  Am  *  Lc^Az 
verified.   From  the  plot,  "TL  441  transfer  characteristics, 
the  log  amplifier  has  1  1/2  decades  of  dynamic  input  range, 
from   20  -  600  mV. 

iii.  Antilog  Configuration 
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Antilog  Configuration  1 


!i 

V_o 

Vo 

+  .028 

-  -OCil 

-  ,o<r^ 

o 

-   .CsXCs 

-  ,  c-z<=» 

,<3«Z 

-,c^ 

-  ,  OJ  V 

.  jx-S 

—  ,i.c*a 

-  ,  ©<*"2 

« *~  s^ 

—  .2CsC» 

,  a'ZV 

,Q<5 

—  .3<^ 

/Ofe^ 

,  uli 

— ,-3CiO 

,i."7L 

.l^H 

-  .<£» 

^S^ 

( 

-  .  £p£i 

^S<* 

SWv^ 

-  ."7<U. 

/z=\^ 

-  irzo 


/ 


^ 


,-svo  .^  ■*-"** 


-88- 


The  antilog  data  is  graphed  in  *'TL  441  antilog 
characteristics"  for  configuration  1.   Configuration  2  will 
expand  the  antilog  dynamic  range  somewhat  above  the  1  1/2 
decades  available  for  configuration  1. 
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F. 


TEST  10 


i.   Purpose 

To  incorporate  the  4  41  log  amplifier  and  the 
quadrant  multiplier  in  the  variable  exponent  circuit, 
general  wiring  is  sketched  below: 


The 


iw^r 
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cc 


-WY 


ii.  Steady  State  Performance 

The  circuit  demonstrated  the  ability  to  track  a 
desired  variable  exponent.   It  was  of  primary  interest 
to  determine  if  the  thermal  properties  were  improved  over 
the  simple  log-fed  multiplier  circuit. 
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The  following  data  was  generated: 
Variable  Exponent  Test 

A..A1  1..AT2. 


2  5  Samples 
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This  represents  a  50%  reduction  in  the  drift  over 
the  simple  log  configuration.   Since  the  441  log  amplifier 
is  fairly  cheap,  it  is  substituted  for  the  simple  log 
configuration . 
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APPENDIX   IV 


DYSYS  SIMULATION 
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APPENDIX   V 


DESIGN  DRAWINGS  AND  SKETCHES 
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APPENDIX   VI 


SPECIFICATIONS 


A. 


ELECTRONIC  COMPONENTS 

1.  Operational  Amplifiers 

a.  Fairchild  a\k   741 

Maximum  ratings 

supply  voltage  +  18v 

Power  dissipation  500  MW 

Supply  current  2 . 8  MA 

Typical  input  offset 

voltage  1  mv 

b.  Fairchild  A<A   32  4  Quad  OP  amp 

Maximum  ratings 
supply  voltage 
single  polarity 

Input  voltage  range 

Power  dissipation 

Supply  current 

Typical  input  offset 
voltage 

2.  EXAR  2208  Four  Quadrant  Multiplier 

Maximum  ratings 

supply  voltage  +  18v 

Power  dissipation  750  MW 

Typical  supply  current  5  MA 

Typical  output  nonlinearity       .  5% 


32v 

-0.3v  to  32v 

900  MW 

2  MA 

2  mv 


+  8v 

6v 

500  MW 

40  Mv 

18.5  MA 

8.  5  MA 
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Texas  Instruments  TL  441  Logarithmic 
Amplifier 

Maximum  ratings 
supply  voltage 

Input  voltage 

Power  dissipation 

Typical  differential 
output  offset  voltage 

Typical  supply  current  vcc 

vcc 


4.   Trimming  Potentiometers 

All  applications:   Spectrol  type  64 
ceramic  potentiometers,  l-100Kft 


Resistors 

All  applications:   Corning  2% 
1/4  watt  resistors 


6.   Power  Supply 

Analog  Devices  model  975  dual  polarity 
power  supply 

Output  voltage  +  15v 

Maximum  rating  500  MA 

Typical  output  error  1% 


B.        PERIPHERALS 


1.   Breadboard  Panels 

Continential  Specialties  model  QT  59 
6.5"  panels 
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2 .   Volt  Meters 

Beede  Instrument  3"  panel  meters, 
model  3-03-8 

o 
250  movement 

200fi  internal  resistance 


3.   Joystick 

2  -  axis  movement; two,  lOOKfi  rotary 
potentiometers  on  each  axis. 


4.   Chassis 


.062"  clear  anodized  aluminum  assembled 
with  pop-rivet  construction. 
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